Magnetization dynamics in dysprosium orthoferrites via inverse Faraday effect 
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The ultrafast non-thermal control of magnetization has recently become feasible in canted 
antiferromagnets through photomagnetic instantaneous pulses [A.V. Kimel et al, Nature 435, 655 
(2005)]. In this experiment circularly polarized femtosecond laser pulses set up a strong magnetic 
field along the wave vector of the radiation through the inverse Faraday effect, thereby exciting 
non-thermally the spin dynamics of dysprosium orthoferrites. A theoretical study is performed 
by using a model for orthoferrites based on a general form of free energy whose parameters are 
extracted from experimental measurements. The magnetization dynamics is described by solving 
coupled sublattice Landau-Lifshitz-Gilbert equations whose damping term is associated with the 
scattering rate due to magnon-magnon interaction. Due to the inverse Faraday effect and the 
non-thermal excitation, the effect of the laser is simulated by magnetic field Gaussian pulses 
with temporal width of the order of hundred femtoseconds. When the field is along the z-axis, a 
single resonance mode of the magnetization is excited. The amplitude of the magnetization and 
out-of-phase behavior of the oscillations for fields in z and -z directions are in good agreement with 
the cited experiment. The analysis of the effect of the temperature shows that magnon-magnon 
scattering mechanism affects the decay of the oscillations on the picosecond scale. Finally, when 
the field pulse is along the x-axis, another mode is excited, as observed in experiments. In this 
case the comparison between theoretical and experimental results shows some discrepancies whose 
origin is related to the role played by anisotropies in orthoferrites. 

PACS numbers: 78.20.Ls, 78.47.+p, 75.30. Ds 
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I. INTRODUCTION 



In the last years the need of enhancing the speed of 
modern spin-electronic and magneto-optic devices, and 
of further developing the magnetic storage technology, 
has stimulated many studies aimed at achieving a fun- 
damental understanding of the mechanisms of magneti- 
zation dynamics and switching. In addition to ultrafast 
magnetic field pulses which require complex devices for 
their generation, the spin dynamics has been induced by 
ultrafast optical laser pulses. Recent experiments have 
shown that significant demagnetization of magnetic com- 
pounds can be measured on the scale of a few hundred 
femtoseconds. 1,2 ' 3,4 Typically the light is absorbed by the 
material, giving rise to a rapid increase of temperature 
responsible for the change of magnetization^Sii and spin 
reorientation^ The question concerning the exact speed 
of the initial sub-picosecond magnetization breakdown is 
still subject of debate and partly relates to the question 
of how to interpret magneto-optical experiments4i2iifl 
Furthermore, the cooling time could limit the repeti- 
tion frequency whose value is fundamental for actual 
applications iii 

Recently, non-thermal ultrafast optical control of mag- 
netization has been achieved in canted antiferromagnet 
samples of dysprosium orthoferrites by using circularly 
polarized femtosecond pulses. 12 Via the inverse Faraday 
effect, the light excitation acts on the spins of the system 
as a magnetic field pulse directed along the wave vec- 
tor of the radiation and proportional to its intensity^ 1 ^ 
The inverse Faraday effect does not rely on absorption 
and has its fingerprint in the fact that the helicity of the 
pump controls the sign of the photo-induced magneti- 



zation. This effect plays a role also in the femtosecond 
photomagnetic switching of spins in ferromagnetic garnet 
films J£ 

Since the manipulation of spins by means of circularly 
polarized laser pulses represents an advance in the field 
of ultrafast magnetization dynamics, we analyze thor- 
oughly the experimental work by Kimel et alH and dis- 
cuss its pecularities. In this experiment the difference 
between the Faraday rotations induced by right- and left- 
handed polarized pulses has been studied in the tem- 
perature range between 20 K and 175 K. A characteris- 
tic spin-wave mode, called quasi-antiferro mode, is ex- 
cited by the light pulse along the z-direction. When 
the sample is heated, the frequency of the mode oscil- 
lations increases and the amplitude decreases from the 
low-temperature maximum value of the order of Ms/16, 
where Ms is the saturation magnetization. We point out 
that there is no theoretical explanation concerning the 
magnitude and the temperature behavior of the oscilla- 
tions. Furthermore, some aspects of the cited experiment 
deserve attention. Actually, in the so called T4 phase, 
stable at temperatures higher than 50 K, the oscillations 
have temperature-dependent frequencies in full agree- 
ment with those measured by Raman experiments^ 1 ^ 
On the other hand, at temperatures below 50 K, the fre- 
quency of the photoinduced magnetization stays constant 
in contrast with the results of the Raman spectra and 
other measurements which signal the discontinuous tran- 
sition to another phase, called T\. Therefore, the experi- 
mental data obtained by laser pulses reflect the excitation 
of resonance modes characteristic of the T4 phase even at 
very low temperatures. Finally, the field pulse is also di- 
rected along the x-axis: another spin-wave mode, called 
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quasi-ferro mode, is excited and is characterized by an 
extraordinarily, not well understood, small amplitude. 

In order to obtain a deeper understanding of the non- 
thermal control of magnetization, in this paper we per- 
form simulations related to the experiment by Kimel et 
alw^k We have studied the magnetization dynamics em- 
ploying a model for orthoferrites that was previously pro- 
posed for the analysis of resonance and high-frequency 
susceptibility^ The parameters of the free energy, such 
as the symmetric and antisymmetric exchange, and the 
anisotropy constants, are determined by using the exper- 
imental Raman spectra of Ref.17. The dynamical behav- 
ior is described by solving two coupled sublattice non- 
linear Landau-Lifshitz-Gilbert equations through a fifth- 
order Runge-Kutta algorithm. The damping term in the 
dynamical equations is related to magnon-magnon inter- 
action and its temperature behavior is provided by a cal- 
culation of the scattering rate in orthoferrites. 22 Exploit- 
ing the inverse Faraday effect, we have analyzed the ef- 
fect of Gaussian magnetic field pulses whose time width 
is of the order of hundred femtoseconds. Since, in the 
regime considered in the experiments, the effective mag- 
netic fields are not large if compared with exchange fields, 
the solution of the linearized system represents a reason- 
able approximation to the numerical results. Therefore 
we have studied the dynamics within the linear solution 
after the excitation by a pulse shaped as a delta function, 
since the magnetic field pulse takes place on a time scale 
shorter than the period of the resonance modes. 

One result of this work is that, in the T4 phase, the 
quasi-antiferro mode of the magnetization is excited by a 
field pulse along the z-axis and the oscillations have am- 
plitudes in agreement with experimental results. More- 
over, the oscillations induced by pulses directed along the 
z and -z axis show the characteristic out-of-phase behav- 
ior. We point out that, even for the dynamics, the ratio 
between the antisymmetric and symmetric exchange en- 
ergies is important. Furthermore, we stress that in the 
canted antiferromagnets, such as rare-earth orthoferrites, 
the largest amplitudes of the oscillation are not obtained 
for the ferromagnetic sum vector of the sublattice mag- 
netizations but for the antiferromagnetic difference vec- 
tor. The behavior of the magnetization has been ana- 
lyzed in the temperature range between 20 K and 175 
K. The damping process based on magnon-magnon scat- 
tering describes the decay of the oscillations with results 
consistent with the experiment. 

The case of the field pulse directed along the x-axis has 
been analyzed in the T4 phase. The magnetization along 
the x-axis oscillates with the frequency of the quasi-ferro 
mode as found in the experiment. However, the calcu- 
lated and the experimental amplitudes of the oscillation 
are different if the magnetic field pulse along the x-axis 
has the same intensity of that along the z-axis. Therefore, 
in the comparison with experimental data, we discuss the 
role of anisotropics between z and x directions as a source 
of discrepancy between theory and experiment. 

Finally, we consider the actual stable phase in equi- 



librium at low temperatures, the T\ phase, and its reso- 
nance modes. We point out that the difference in energy 
between the Fx and T4 can be very small, so that even 
a small laser-heating effect could be responsible for the 
stabilization of the T4 phase at very low temperatures on 
a picosecond time scale. 

The outline of this paper is as follows. In the next 
section we discuss the numerical approach for the sys- 
tem and the analytic solution of the linearized equations 
for excitation by a delta function shaped magnetic field 
pulse. Section III provides the numerical and analytical 
results: in the first and second subsection the excitation 
due to the pulses along the z-axis and x-axis, respectively, 
is analyzed when the system at equilibrium is in the T4 
phase. In the final subsection the effect on the dynamics 
of a Ti phase stable at low temperatures is considered. 
Section IV provides a summary. 



II. FREE ENERGY AND DYNAMICAL 
EQUATIONS 

Rare-earth orthoferrites are represented by the formula 
ReFeO?,, where Re stands for rare-earth. They have a 
perovskite-type structure with slight deformation from 
cubic to orthorombic. In many of them the spins of iron 
ions are antiferromagnetically aligned through a strong 
super-exchange interaction with a Necl transition tem- 
perature of about 700K. Moreover, promoted by the or- 
thorombic deformation, an antisymmetric exchange in- 
teraction acts between iron spins, resulting in spin canted 
magnetism with a feeble saturation moment. In the case 
of dysprosium orthoferrites, the temperature dependence 
of the ferromagnetic moment is characterized by a steep 
rise around 50 K in coincidence with the stabilization 
of the T4 phased Actually these compounds are also 
known for their spin reorientation properties: continu- 
ous rotational-type (ferromagnetism present in the low- 
temperature phase) and, only for dysprosium, abrupt- 
type from Ti (antiferromagnetic) to T4 with increasing 
temperature. 

We use the free energy and the dynamical equations 
proposed in a previous work and focus on the behavior 
of the magnetization in the T4 phased The static and 
dynamical behavior is studied within the two-sublatticc 
model for iron spins that takes into account the two ac- 
tive resonance modes. This represents an excellent ap- 
proximation since the remaining resonance modes are al- 
most inactive, and their interaction with active modes is 
negligible. In Fig.l we show a schematic representation 
of equilibrium positions of two-sublattice magnetization 
vectors for dysprosium orthoferrites^ 

The normalized free energy V = F/Mq, with Mo the 
modulus of the sublattice magnetization, is composed of 
a part V exc due to exchange interactions and a part V an i 
due to the anisotropy: 

V = V exc + V am . (1) 
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FIG. 1: The equilibrium positions of sublattice magnetization 
vectors for the high-temperature Fi (ferromagnetic) and low- 
temperature Fi (antiferromagnetic) phases. Due to the ra- 
tio between the antisymmetric D and symmetric E exchange 
fields, the angle /3o is actually very small. 



The free energy is expanded as a power series in the mag- 
netization components and, in order to properly describe 
the T4 phase, only quadratic terms are sufficient. The 
exchange energy is written as the sum of a scalar and a 
pseudo- vector part 



ERi ■ R 2 + D(X X Z 2 - X 2 Z X ) 



(2) 



where E and D are the symmetric and antisymmetric ex- 
change fields, respectively, -Ri = Mi/M = (X\,Yi,Z\), 
and R 2 = M 2 /M = (X 2 ,Y 2 ,Z 2 ). The anisotropic en- 
ergy is 



V ar , 



A XX (X\ + Xl) - A zz {Zl + Z\). 



(3) 



The exchange field E is quite large compared to the other 
terms. It is related to the exchange spin-spin interaction 
J via the equation E = l2JS/gfiB, with S = 5/2 the spin 
of iron ions, g — 2, and \ib the Bohr magneton. Since J 
is of the order of 20 cm -1 , E is approximately 6.4 x 10 6 
Oeiii The field D is related to the antisymmetrical ex- 
change energy d through the relation D = 6dS/g[iB and 
is of the order of 1.4 x 10 5 Oe (d is about 0.88 cm -1 ). 
The anisotropy constants A xx and A zz depend on tem- 
perature and are of the order of hundreds of Oersted. 

Using the energy given in Eq.Q), the equilibrium po- 
sition can be derived for the T4 phase. This is shown in 
Fig. 1, where the small canting angle (3q is determined 
by the equation 



tan(2/3 Q ) = 



D 



A, 



D 

— = 0.022. 

E 



(4) 



Hence the magnetization Mg — \M X + M 2 \ = 
2Mosin(/3o) ~ 0.022A/o is two orders of magnitude less 
than M Q . 



The equilibrium position of M\ and M 2 represents 
a stationary solution of the nonlinear Landau-Lifshitz- 
Gilbert equations 



1 dR 1 

1 dR 2 
j~dT 



= A (S(t) - ViV) + <xRi A 
= -R 2 A [S{t) - V 2 V 



R 2 A 



au 2 



dR 1 
~dT 



dR 2 
dt 



(5) 



(6) 



where 7 = 17.6 MHz/Oe is the gyroscopic ratio, Vi and 
V2 are gradients with respect to R\ and R 2 , respectively, 
and V is the energy in Eq.Q. Clearly the dynamical 
equations satisfy the following constraints: X\ + Y± + 
Z\ = 1 and Xl + Y 2 2 + Zf = 1. 

The quantity H(t) in Eas. (|5l6J) is the time-dependent 
magnetic field simulating the effect of laser pulses due to 
the inverse Faraday effect. In the experiment by Kimel et 
alm& the magnetic field pulse is of the order of fractions 
of Tesla with time width of hundred femtoseconds. In the 
numerical simulations we consider a pulse directed along 
the propagation direction of the light with a Gaussian 
shape 



H{t) 



, F 



exp [-{t/T p ) 2 ] , 



(7) 



where k defines the direction of the light wave- vector and 
t p indicates approximately the duration of the pulsed 

Finally, in Eas. (|5l6(l a is the Gilbert constant. Actu- 
ally, a takes into account the damping of the oscillations 
due to the magnon-magnon scattering and to the interac- 
tion of magnons with dysprosium spins and phonons. We 
notice that the scattering via dysprosium spins should 
be larger at very low temperatures where dysprosium 
ions tend to order. Moreover, the phonon-magnon scat- 
tering should be effective only on the nanosecond time 
scale^*2i It is the magnon-magnon interaction that pro- 
vides the preminent source of scattering on the picosec- 
ond time scale. 



A. Solution of linearized system 

In this subsection we consider the solution determined 
by linearizing Eqs.(0EJl and study the excitation of this 
linear system due to a magnetic field pulse shaped as a 
delta function. This is reasonable since, in the regime 
considered in the experiments, the effective magnetic 
fields are not large if compared with exchange fields, and 
the temporal width of the pulse is much smaller than the 
periods of the resonance modes. 

In order to take into account small deviations from 
the equilibrium, the standard approach is to con- 
sider two separate coordinate systems, (Si,Ti,Y±) and 
(S 2 ,T 2 ,Y 2 ), which describe the dynamics of M\ and M 2 , 
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respectively. 18 The variables Si and S 2 are chosen in or- 
der to coincide with the equilibrium positions of M\ and 
M 2 , respectively, so that 



and 



S 1 ! = sin(/? )Zi +cos(/? )Xi, 
Ti = -cos(/3 )^i+sin(/3 )X 1 , 

52 = sin{f3 a )Z 2 - cos(/3 )^2, 
Ti = cos(/3 )^2+sin(/3 )X 2 . 



(8) 
(9) 



(10) 
(11) 



By linearizing the system, we obtain the frequencies of 
two modes, the quasi-antiferro loafm and the quasi-ferro 
lofm modes, involving cooperative motions of spins of 
the two sublattices. The energy of the quasi-antiferro 
and quasi-ferro modes is of the order of several cm -1 . 
The first mode is characterized by the frequency 



W AFM 
7 2 



AEA XX + 4A XX (A XX - A zz ) + D 2 



(12) 



The dynamic in this mode shows the following behav- 
ior: AX x {t) = -AX 2 (t), AYi (t) = -AY 2 (t) and 
AZt(t) = AZ 2 (t), where AWi(i) = W l {t) - W- q , with 
W = X, Y, Z, i = 1,2, and W^ 9 corresponding to the 
equilibrium position shown in Fig.l. Thus, the only fer- 
romagnetic component different from zero is the magne- 
tization along the z-axis with respect to the equilibrium, 
AM Z = MqAMz, where AM Z is defined by 



AM Z = AZdt) + AZ 2 (t) = Z x [t) 



z 2 {t)-{zl q + z e 2 q ). 

(13) 

Since the net spin does not reflect the motions of the sub- 
lattice spins, this mode is called quasi-antiferromagnetic. 
The second mode, with frequency 



FM 
7 2 



AE{A XX - A ZZ ) + AA XX {A XX - A zz ), (14) 



is characterized by AX^t) = AX 2 (t), AFi(i) = AY 2 (t) 
and AZi(t) = —AZ 2 (t). Since the net spin executes the 
same rocking behavior as the sublattice spins, this mode 
is called quasi-ferromagnetic. 

Eas. (|12ll4f> relate the mode frequencies to the model 
parameters. The exchange fields are assumed constant 
in temperature at the values given above, since they rep- 
resent the highest energy scales. Using Eas. (|12ll4l) and 
the experimentally measured Raman spectra of Ref.17, 
we derive the anisotropy parameters. The quasi-antiferro 
mode increases linearly as a function of temperature in 
the T 4 phase: from 150 GHz at 50 K to 450 GHz at about 
200 K, while the quasi-ferro mode stays constant at about 
375 GHz. Therefore the anisotropic terms change upon 
heating the sample: A xx varies from about -640 Oe at 
T=50 K to about 200 Oe at T=200 K, A zz from -1540 Oe 
to -700 Oe in the same temperature range, with approxi- 
mately fixed A xx — A zz = 900 Oe. The knowledge of the 
frequencies is important also to determine the value of the 



damping constant in the dynamical equations l1">l(il) . In- 
deed, if ojq is the temperature-dependent frequency of one 
of the modes, the damping constant a could be related 
to the damping rate 1/tq by the relation 1/tq = luqch. At 
the picosecond scale, the dominant spin-wave damping 
is due to four-magnon scattering. By using many-body 
perturbation theory, the rate for spin-wave at zero wave- 
vector is estimated to be l/r ~ 2.66 • 10 4 T 2 s _1 , with T 
temperature in units of Kelvin degrees^ The scattering 
rate grows quadratically in temperature, while the fre- 
quency of the quasi-antiferro mode is an increasing linear 
function of the temperature^ The quantity a gets larger 
with temperature since the most important contribution 
is given by 1/tq. Actually, the values of a corresponding 
to the quasi-antiferro mode range from about 0.4 • 10~ 4 
at T=50 K to 3 • 10~ 4 at T=200 K. The smallness of a 
implies that the oscillations of the magnetization are not 
strongly damped. 

It is useful to consider a magnetic field pulse shaped 
as a delta function 



H(t) = kF S{t). 



(15) 



Since the pulse is instantaneous, it provides an initial 
condition to the dynamics described by the linearized 
equations of motions. We have analyzed two cases: field 
along z- and x-axis, since these are prominent for the 
experiment that we want to discuss. 

For the field along the z-axis, starting from the equi- 
librium position at t = 0~, we find at t = 0^ 



AXj^ 
AYx(0 4 
AZj^ 



= -AX 2 (0+) = cos^oMcos^Fo) - 1](16) 
- -Ay 2 (0+) - cos(/3 ) sin( 7 F ) (17) 
= AZ 2 {0+) = 0. (18) 



Therefore this pulse will excite only the quasi-antiferro 
mode. 

The field along the x-axis yields at t = + 



AXi(0+) 
AF!(0+) 
AZ!(0+) 



= AA 2 (0+) = (19) 
= AY 2 (0 + ) = -sin(/3 )sin( 7 Fo) (20) 
= -AZ 2 (0+) = sin(/? )[cos( 7 F ) - 1]<21) 



Thus the system will be excited in the quasi-ferro mode. 

As the experimental field magnitude is lower than the 
exchange fields, we study the subsequent dynamical evo- 
lution within the linear approximation. 18 Within a relax- 
ation time approximation, we include the damping time r 
describing the decay towards the equilibrium position. 23 

From the solution of the linear system in the case of a 
pulse along the z-axis, we get the evolution of the quasi- 
antiferro mode. The only ferromagnetic component, de- 
fined in Eq.{T3J, is 

AMz = e -'/r£2!^l sin ( 2/ 3 ) sin 2 (^) co S (w AFM t) 

+2e"* /r R cos 2 (j3 ) sin(7F ) Bxa.(<jj AFM t) 
+2e-*/ r sin/3 cos 2 (/3 ) [cos( 7 F ) - 1] . (22) 
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We notice in Eq. ill'L'l) that the even terms in the magnetic 
field are at least two orders of magnitude less than the 
odd term. 

Along the other two directions, there are the antifer- 
romagnetic component along x-axis (with respect to the 
equilibrium), A Ax — Mo A Ax, with AAx defined as 

AA X = AX 1 (t)-AX 2 (t) = X 1 {t)-X 2 {t)-(X e 1 q -X e 2 q ), 

(23) 

and the antifcrromagnetic component along y-axis (with 
respect to the equilibrium), AAy = Mo AAy, with AAy 

AAy = AY 1 (t) - AY 2 (t) = Yx(i) - Y 2 (t) - (Y* 9 - Y 2 eq ). 

(24) 

Using again the linearized form of Eas. i|5l6[) . we obtain 



AA 



v : -e- t/T sin(2/3 ) sin(/3 ) sin 2 (^) cos(u AFM t) 

-e _t/r i?sin(2/3 ) sin(7-F ) s\n{uj AFM t) 
+2e- f / T cos/3 cos 2 (/3 ) [cos( 7 F ) - 1] , (25) 



and 
AA Y 



-e-V^I sin(2/? ) sin 2 (^) sia(cj AFM t) 
+2e' t/T cos(/3 ) sin(7F ) cos(ui AFM t). (26) 



The components specified in Eas. l|22l25l26|) are deter- 
mined by the mode frequency loafm, the angle (3q, and 
the quantity R which depends on the model parameters 
in the form 



R = 



2A XX + D 2 /2E 
2E 



(27) 



At intermediate temperatures (A xx very small) one gets 
R-D/2E. 

When the pulse is directed along the x-direction, the 
quasi-ferro mode is excited. The components of the mag- 
netization different from zero are those along x- and y- 
axis. In particular we consider the magnetization along 
x-axis (with respect to equilibrium), AAIx = MqAMx, 
where AMx is defined as 



AM X — AXi{t) + AX 2 (t). 



(28) 



From the linearized dynamical equations this quantity is 
calculated as 



AM 



x 



III. 



-t/r 



AE 



A x 



A z 



■ sin 2 (/3 ) sin(7F ) sa\(uj FM t). 

(29) 



NUMERICAL AND ANALYTICAL 
RESULTS 



Using the procedure outlined in the previous section, 
all the parameters appearing in the dynamical equa- 
tions lloltiH can be determined. These nonlinear equa- 
tions are numerically integrated through a fifth-order 



Runge-Kutta algorithm^ The experimental laser pulses 
along the z-axis are estimated to be equivalent to mag- 
netic fields with an amplitude of 0.3 T and a full-width 
at half- maximum of about 200 fsii 2 " Therefore, we con- 
sider Gaussian fields of the form given in Eq.Q, with 
Fq = ^TTTpAmp(H), t p = 200 fs and Amp(H) typically 
3000 Oe. This value of the effective field is small when 
compared with the exchange fields, so that effects due 
to the nonlinearity of the equations are negligible. Since 
the time width of the pulse is small on the scale of the 
mode periods, the results are not strongly dependent on 
T p . The linear solution with a delta pulse represents a 
good approximation, provides the right orders of mag- 
nitude but does not reproduce exactly the numerical re- 
sults. Therefore, in all the following figures we plot the 
results obtained by the numerical integration. 

In the first and second subsections we discuss the ef- 
fects of a Gaussian pulse along z and x directions, respec- 
tively, when the system at equilibrium is in the I^ phase. 
In the final one we briefly analyze the resonance modes 
and magnetization dynamics if the Ti phase is stable at 
low temperatures. 



A. Pulse along z-direction 

As illustrated in Ref. 12, pulses along the z and - 
z directions excite the quasi-antiferro mode. Using the 
numerical procedure discussed in the preceding section, 
we have calculated the oscillating behavior of the ferro- 
magnetic vector along z focusing on AMz at T = 95 
K. As shown in Fig. 2(a), the magnetization has a sine- 
like behavior with weak damping on a picosecond scale. 
In addition, opposite fields give rise to out-of-phase os- 
cillations in agreement with the experimental data. We 
notice that after 50 picoseconds the amplitude is about 
half of the initial one, a value that is compatible with re- 
sults reported in Fig.l of Ref.12. This suggests that the 
damping term associated to magnon-magnon scattering 
provides a reasonable description of the reduction of the 
amplitude as a function of temperature. 

Comparing the experimental data shown in Fig.l of 
Ref.12 with the theoretical results shown in Fig. 2(a) of 
this paper, we point out similarities and differences. In- 
deed, on the scale of the initial pulse width, there is a 
strong enhancement of the experimental Faraday rota- 
tion, probably due to the interference between pump and 
probe pulses. After this transient the oscillations induced 
by pulses with opposite helicities show similar amplitudes 
in time. In contrast with theoretical results, the equilib- 
rium point of the Faraday rotations shows a decay and 
its behavior is different for left- and right-handed po- 
larized pulses. Only later the oscillations tend toward 
a common equilibrium point that is different, however, 
from that before the pulse excitation. If the shift of 
the equilibrium positions is caused by a change of the 
anisotropy constants due to an intrinsic photomagnetic 
effect, then the amplitude of the oscillations should be 
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different for the two helicities, as reported in a recent 
experiment. 15 Since this is not the case in dysprosium 
orthoferrites, the change in the orientation equilibrium 
is more probably associated with a small but unavoid- 
able laser-heating effect. Actually, within the model pro- 
posed in this paper, a shift of the equilibrium configura- 
tion could in principle be caused by the magnetic pulse 
during the excitation^ This shift is negligible, however, 
since the amplitude of the effective magnetic field is much 
smaller than the antisymmetrical exchange field D. Only 
by enhancing the amplitude of the pulse to values com- 
parable to D, the shift becomes important, and, upon 
further increase, even a precessional switching can occur. 
These amplitudes imply very high light intensities that 
are probably above the damage threshold of the samples. 

In Fig. 2(b) we focus on the related quantity d(AMz), 
defined as 

d(AM z ) = AM Z (H+) - AM Z (H-), (30) 

with H+ and H— indicating positive and negative am- 
plitudes, respectively. As reported in the experimental 
workfi^ this quantity is less dependent on initial effects. 
Just after the transient induced by the pulse, this quan- 
tity has an amplitude of the order of Ms/16, in close 
agreement with the calculated data reported in Fig. 2(b). 
It is worthwhile understanding the order of magnitude 
of this amplitude by exploiting the result of the previous 
section for the delta pulse. From Ea. H22|) we get for short 
times (t <C T2) 

AM Z ~ 2A/ i?sin( 7 Fo) rin(u AFM t), (31) 

with R given in Ea. (|27|l . The amplitude is determined 
by the term R depending on the model parameters and 
by sin(7i 7 o) related to the pulse intensity. Using the ex- 
perimental data, the impulse Fq is estimated to be of the 
order of 3000 x 400 Oe-fs. This implies that 7.F0 is of 
the same order as the angle (3q responsible for the canted 
antifcrromagnetism: sin(7_fo) ~ 7-F0 ~ 2/3q. Thus, we 
have 

AMz/Ms ^ 2Rsm{u AFM t). (32) 

If we consider R ~ D/2E, we obtain the order of mag- 
nitude lQd(AM z )/M s - 16 x 2 x 0.02 = 0.64 not so far 
from the experiment. Hence the ratio D/E is fundamen- 
tal not only for the equilibrium configuration, but also 
for the magnetization dynamics in the canted antiferro- 
magnet. 

Within the theoretical approach, the dynamics of the 
antifcrromagnetic vectors can be easily obtained, but 
these quantities have not been measured in the exper- 
imental work^2 Due to the fact that the equilibrium po- 
sition corresponds to a maximum of Ax — Mq(X\ — X2), 
the temporal evolution of A Ax = Mo A Ax, with A Ax 
defined in Ea. (|23|l . is characterized by a very small ampli- 
tude. From Ea. (|25|l one can deduce that the amplitude 
of AAx should be at least an order of magnitude smaller 
than that of AM Z . 



0.3 




FIG. 2: (a) The variation AMz of the magnetization along 
z-axis with respect to the equilibrium value in units of Ms/16 
as a function of time for two opposite amplitudes, (b) The 
difference d(AMz) in units of Ms/16 between the magne- 
tizations obtained upon excitation by opposite fields along 
z-direction as function of time. The dashed line includes the 
effect of the anisotropy that will be introduced in the next 
subsection. 

In Fig. 3 we show the results of the numerical calcu- 
lation for the quantity AAy — Mo AAy, where AAy is 
defined in Eq. I|24|) . We point out that now the amplitude 
is an order of magnitude larger than that of AMz an d 
the response is more sensitive to the pulse. Actually this 
is due to the fact that the dynamics of the components 
along the y-axis is strongly influenced by the highest en- 
ergy scale From Eq. Q26J1 we obtain a rough estimate 
for short times (t < T2) 

AA Y 2M sin(7F ) cos(u} AFM t). (33) 

Therefore the amplitude is 2sin( 7 Fo)Mo ~ 2(3qMq ~ 
0.022M = Ms. Finally, we notice that, after 50 pi- 
coseconds, the damping acts in the same way as for the 
ferromagnetic vector along the z-direction, causing a re- 
duction of about half of the initial amplitude. 

In order to make contact to the experimental measure- 
ments, it is interesting to analyze the behavior of the 
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FIG. 3: The variation of the antiferromagnetic vector along 
y-axis AAy with respect to Mo as a function of time for two 
opposite amplitudes. 



magnetization dynamics at different temperatures. As 
reported in Fig. 3 of the paper by Kimel et al.^ upon 
heating the sample, the frequency of the oscillations in- 
creases and the amplitude decreases. The oscillations 
have temperature-dependent frequencies that are very 
close to those measured by Raman experiments at tem- 
peratures higher than 50 K. In Fig. 4 of this paper the 
calculated d(AMz) is shown. We find agreement with 
experimental decay of amplitudes and behavior of the 
frequency upon increasing the temperature for times af- 
ter the transient. Within the theoretical approach, the 
temperature dependence of the anisotropy constants A xx 
and A zz affects the frequency of the modes, but is not 
so important in the change of the amplitude after the 
laser transient. The decrease of the magnetization seen 
at a fixed time (of the order of tens of picoseconds) for 
different temperatures is dominated by the increase of 
the damping constant a in temperature. Due to the 
agreement with experiment, the estimate of the damp- 
ing constant derived on the basis of the magnon-magnon 
scattering is reliable. 

The results of Fig. 4 show the magnetization normal- 
ized by the sublattice magnetization modulus Mo. We 
point out that, in principle, Mq can vary as function of 
temperature. However, the inclusion of the temperature 
dependence of Mq — M z /(2 sin(/?o)) is not easy, since 
only the magnetization M z is typically experimentally 
available. In any case, even if the angle [3q is assumed 
fixed, the change of M z is small in the considered tem- 
perature ranged* 

Before closing this section, we focus on the dependence 
of the magnetic response on the amplitude of the pulse. 




50 100 150 200 



t(ps) 

FIG. 4: The difference between the magnetizations along z 
for opposite fields as a function of temperature. Every new 
curve is shifted from the previous one along the vertical axis 
over 1. 

As inferred by Eq. i|31|) . for the effective magnetic pulses 
used in the experiments, the ferromagnetic vector is pro- 
portional to Fq, *- e -i to the intensity of the light pulse 
used in the inverse Faraday effect. In Fig. 5 we plot the 
numerical results for the amplitude of the magnetization 
at t=40 ps, i.e., after about 10 periods as in the exper- 
imental data^ We find the expected linear behavior as 
a function of the amplitude of the pulse. The value ob- 
tained at 5 T for T — 60 K is not far from unity. Thus, 
the saturation of the magnetization along z is reached for 
this high intensity, in agreement with the extrapolation 
of the experimental data in the inset of Fig. 2 of Ref.12. 

Summarizing, in this subsection we have focused on the 
magnetization dynamics in the T4 phase when the field 
pulse is directed along the z-axis. We have analyzed the 
amplitude and the decay of the oscillations as a function 
of temperature finding agreement with experimental data 
for times longer than the initial transient. Actually, on 
the time scale of the initial pulse width, discrepancies 
with experiments appear. These could be due to a small 
but unavoidable laser-heating effect that is not included 
in the present theoretical approach. 



B. Pulse along x-direction: role of anisotropy 

In this subsection we analyze the effects of the excita- 
tion due to the pulse along the x-axis. 

The fields along the x-direction excite another mode. 
As shown in Ea. (|29|) . AMx oscillates with the frequency 
of the quasi- ferro mode in agreement with experiment. 
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Amp(H) (10* Oe) 

FIG. 5: Amplitude of the spin oscillations as a function of the 
amplitude of the magnetic field pulse for different tempera- 
tures. 



On the other hand, the magnitude of the magnetic re- 
sponse along x is of the same order as that obtained 
along the z-axis with a pulse of equal amplitude. Even 
if the torque exerted by the field along x is tiny, the 
subsequent temporal evolution, in particular, that of the 
variables T\ and T2, is able to give a non-negligible am- 
plitude to the response along x. This is in discrepancy 
with the experimental results which show that the Fara- 
day rotation along x is at least an order of magnitude 
smaller than that along the z-axisii£ Also, as shown in 
Fig. 3 of Ref.12, it does not show any temperature de- 
pendence following the behavior of the frequency of the 
quasi-ferro mode, in striking contrast with the case of a 
pulse along the z-axis. 

In order to properly compare the experimental and 
theoretical results, we should take into account several 
effects in orthoferritesi2& First of all, the experimental 
estimate of the effective magnetic field pulse due to the 
inverse Faraday effect is available only when the pulse 
is along the z-axis. Moreover, it is well known that 
in rare-earth orthoferrites optical birefringence is not 
negligiblei2L2& Therefore the experimental Faraday ro- 
tations along the z- and x-axis can be different, even if 
the oscillations of the magnetizations are of the same or- 
der. Since the Faraday rotation along the x direction is 
very small in comparison with that along the z-axis, part 
of the effect could be also associated with the anisotropy 
of the magneto-optical susceptibility. Indeed, the effec- 
tive magnetic field generated via the inverse Faraday ef- 
fect can be different for the two orthogonal orientations. 
Finally, there is another source of anisotropy: the field 
along the z-axis can show a renormalization of the cou- 
pling to the system that is different from that along the 
x-axis. 

In order to elucidate this last point, we have included 




10 20 30 40 50 
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FIG. 6: Comparison between amplitudes for excitations with 
fields along z and x axis as function of time at T=60 K. In 
the inset the ratio of these amplitudes as function of the tem- 
perature for times just after the pulse. 



in our model the anisotropy induced by the dysprosium 
ions. As a result, a modification of the coupling of the 
magnetic field to the iron ions takes place. According 
to a simple scheme proposed long ago by Zvezdin and 
Matveev^S the coupling to the field along the z-direction 
is reduced by the factor 1 + rj z ~ 0.74 (see in Fig. 2(b) the 
effect of this reduction). Instead, the reduction 1 + r) x 
of the coupling along the x-axis is temperature depen- 
dent and can be extracted by the measurements of mag- 
netic properties^ With temperature decreasing toward 
50 K, these coupling factors can become a non-negligible 
source of anisotropy. As shown in Fig. 6, the difference 
in the amplitudes for the response along the z- and x-axis 
can become relevant at those temperatures. However, as 
shown in the inset of Fig. 6, the ratio between the am- 
plitudes along x and z axis increases with temperature 
up to an inversion point. Therefore the source of the 
anisotropy associated only to dysprosium ions seems not 
to be effective for the interpretation of the experimental 
data. 



C. Stable phase at low temperatures 

In this last subsection we focus on the Ti phase that 
is stable at low temperatures. 

The measurements reported in Ref.12 show the max- 
imum value of the amplitude of the magnetization be- 
tween 20 and 50 K. Moreover, the frequency of the pho- 
toinduced magnetization is constant in this temperature 
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range at the value close to 50 K that is characteristic of 
the T 4 phase. Therefore, the laser pulse is inducing ex- 
citation modes of the phase that should be unstable at 
low temperatures. 

In order to stabilize the Ti phase, the free energy (JTJ 
has to be supplemented with a quartic term, 29 for in- 
stance, a term that decreases the energy when the mo- 
ments are antiferromagnetically aligned along the y axis 
(see Fig.l). We assume V -► V = V - A {2 J(Y 1 - F 2 ) 4 in 
order to describe the system in simple terms. Upon lin- 
earization around the equilibrium configuration Y^ q = 1 
and Y 2 eq = — 1 > the two modes characteristic of this phase 
can be derived. The first one, with frequency 
, ,2 
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2 



= 2E(32A$ - 2A XX ) - D 2 , 



(34) 



corresponds to a dynamic with AX 4 = —AX 2 and 
AZi = AZ2, while the second mode, with frequency 

^ = 2E(2A zz -32A y 2 J)-D 2 , (35) 

is characterized by AXi — AX2 and AZi = —AZ 2 - 
We follow the procedure of Sec. II: the anisotropy con- 
stants are derived by Eas. (|34l35|l and the experimental 

Raman spectra-! by imposing Ayy — 0. At T=50 we find 
A^Ti) = -1040 Oe that is close to A XX (T 4 ) = -640 Oe 
obtained if the T4 phase is stable at that temperature. 
The values of Ayy required to stabilize the Ti phase are 

small. The free energy corresponding to the equilibrium 

(2) 

position of the Ti phase is Vi = —E — WA yy , while that 
of the T4 phase is at first order in the canting angle (3q 

D 2 



V 4 



-E - D(2(3 ) - 2A 

The condition V\ < V4 implies that 

D 2 A, 



lift} ' -EJ 2^/4.^^ a 

E 



A (2) > 

yy 



WE 



(36) 



(37) 



Taking into account the values of the exchange fields 
given in the previous section and the anisotropy constant 

A xx (Ti), we derive A vy > 60 Oe. Therefore the value of 

(2) 

the constant A yy is consistent with the expansion of the 
free energy being only a fraction of the anisotropy ener- 
gies A xx and A 

zz • 

If the parameter A y 2 J is of the order of hundreds of 
Oersted, the phases T 4 and Ti are close in energy. Due 
to a small light absorption, the laser pulse could af- 
fect the stability of the system by favoring the T 4 phase 
on a short time scale. Therefore the femtosecond pulse 
could induce a reorientational phase transition from Ti 
to T 4 state in analogy with the antiferromagnetic-to- 
ferromagnetic phase transition induced by heating with a 
laser in FeRh films. 30 Finally, we point out that a static 
magnetic field gives rise to a spin reorentational tran- 
sition in orthoferritesi^S Hence the role of the magnetic 
field pulse induced via the inverse Faraday effect could be 
investigated in relation to the perturbation of the phase 
stability. This is left for future investigations. 



IV. SUMMARY 

Stimulated by recent experimental results showing ul- 
trafast non-thermal control of magnetization by instanta- 
neous photomagnetic pulses in dysprosium orthoferrites, 
a theoretical study of magnetization dynamics has been 
presented in this paper. We have employed a general 
form of free energy suitable for dysprosium orthoferrites 
whose parameters are derived from experimental mea- 
surements. We have solved coupled sublattice Landau- 
Lifshitz- Gilbert equations whose damping parameter is 
determined by considering the scattering rate due to 
magnon-magnon interaction. Due to the inverse Fara- 
day effect, the magnetic fields perturbing the equilib- 
rium configuration can be modeled as Gaussian pulses 
with amplitude proportional to the intensity of the light 
pulse and time width of the order of hundred femtosec- 
onds. The non-linear dynamical equations have been in- 
tegrated through an optimized Runge-Kutta algorithm 
and an analytical solution of the linearized system has 
been discussed in the case when the magnetic field pulse 
is assumed to have the shape of a delta function. This 
solution provides the right orders of magnitude allowing 
to interpret the experimental results in simple terms. 

We have found that the quasi-antiferro mode is excited 
by the pulse along the z-axis and the oscillations of the 
magnetization have amplitudes compatible with experi- 
ment. Magnetic fields in opposite directions give rise to 
out-of-phase oscillations showing a behavior in agreement 
with experimental results for times longer than the initial 
transient. We have stressed that the magnetization dy- 
namics is not only strongly influenced by the amplitude 
of the magnetic field pulse, but also by the parameters 
determining the free energy, in particular the ratio be- 
tween the antisymmetric and symmetric exchange ener- 
gies. The temperature dependence of the magnetization 
dynamics has been discussed showing that the proposed 
damping mechanism based on magnon-magnon scatter- 
ing can be effective on the picosecond scale. When the 
field pulse is along the x-axis, the quasi-ferro mode is ex- 
cited, but there are some discrepancies in the comparison 
between theory and data. We point out that the response 
along the x-axis can be strongly influenced in orthofer- 
rites by several effects such as the optical birefringence 
and the anisotropy of the magneto-optical susceptibility. 
Finally, the behavior of the magnetization has been an- 
alyzed in the low-temperature range where, due to an 
unavoidable heating effect, the laser pulse could perturb 
the stability between T\ and T 4 state. 

We notice that the model proposed in this work has ne- 
glected dipolar contributions because they are orders of 
magnitude smaller than exchange and anisotropy terms. 
Moreover, due to the fact that the effective magnetic field 
obtained through the inverse Faraday effect shows spatial 
variations negligible on the microscopic scale, only the 
spin-wave modes at zero wave-vector are excited. Since 
the static magnetization changes slowly in the investi- 
gated temperature range, the presence of spin- waves with 



10 



wave-vectors different from zero should not provide siz- 
able contributions to the dynamic behavior. Therefore, 
the macrospin approximation employed in this paper can 
be considered reliable. 

Finally, we point out that the approach employed for 
dysprosium orthoferrites can be also generalized to de- 
scribe the magnetization dynamics of other rare-earth 
orthoferrites, at least in the phase. The anisotropy 
constants are the only quantities strongly dependent on 
the rare-earth ion, but these do not play a major role in 
affecting statics and dynamics. On the other hand, the 
most important values of the exchange fields are of the 
same order in several rare-earth orthoferrites^ Clearly 
the approach proposed in this paper is suitable for mag- 
netic dielectrics and not for metallic itinerant magnets. 



However, up to now, due to the unavoidable light absorp- 
tion, it has been impossible to ascertain the role played 
by the inverse Faraday effect on the magnetization dy- 
namics of itinerant magnets. 
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